Polydomain structure and its origins in isotropic-genesis nematic elastomers 
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We address the physics of nematic hquid crystaUine elastomers randomly crosslinked in the 
isotropic state. To do this, we construct a phenomenological effective replica Hamiltonian in terms 
of two order-parameter fields: one for the vulcanization, the other for nematic alignment. Using 
a Gaussian variational approach, we analyze both thermal and quenched fluctuations of the local 
nematic order, and find that, even for low temperatures, the macroscopically isotropic polydomain 
state is stabilized by the network heterogeneity. For sufficiently strong disorder and low enough 
temperature, our theory predicts unusual, short-range oscillatory structure in (i.e., anti-alignment 
of) the local nematic order. The present approach, which naturally takes into account the compli- 
ant, thermally fluctuating and heterogeneous features of elastomeric networks, can also be applied 
to other types of randomly crosslinked solids. 

PACS numbers: 64.70.pp, 61.41. +e 



I. INTRODUCTION 

Liquid crystal elastomers [T] are fascinating materials, 
in which the liquid crystalline order is strongly coupled 
to the elasticity of the underlying elastomeric network. 
This gives rise to novel properties, such as soft [H |3] 
and semi-soft elasticity [H [S] , which is of interest in fun- 
damental research and, furthermore, has significant po- 
tential for technological applications. From the perspec- 
tive of statistical physics, these materials, like any other 
randomly crosslinked systems, pose important concep- 
tual challenges, as their physical properties do not de- 
pend just on the conditions under which various physi- 
cal quantities are measured. They also depend on both 
the conditions under which the system is prepared and 
the network heterogeneities necessarily generated via the 
random crosslinking preparation process. A coherent 
theory that takes all these factors into account requires 
two statistical ensembles, i. e. the preparation ensemble 
and the measurement ensemble. This essential point was 
qualitatively understood by Deam and Edwards and by 
de Gennes as early as the 1970's (see, e.g., Refs. |H[7]). 
Its full significance however has not been explored previ- 
ously. 

One such challenge concerns the origins, structure and 
elasticity of the polydomain state in isotropic-genesis 
nematic elastomers (IGNEs; i.e., elastomers that were 
crosslinked in the isotropic state) [8-18 . It has always 
been observed experimentally that IGNEs at low tem- 
perature are characterized by a polydomain (PD) state, 
i.e., a macroscopically isotropic state consisting of well- 
defined domains of randomly oriented nematic order. 
The typical domain-size is of the order of microns. Fur- 
thermore, PD IGNEs exhibit much softer stress-strain 
curves than PD nematic elastomers crosslinked in ne- 
matic state, as demonstrated by recent experiments of 
Urayama's group |16j . 



A heuristic argument, of the type first developed by 
Imry and Ma |f5] for systems at low temperature hav- 
ing quenched random fields with short-range correlations, 
predicts that long-range order involving continuous sym- 
metry breaking is precluded in systems of spatial dimen- 
sions fewer than four. This argument surely applies to 
nematogens in the presence of an immobile random sub- 
strate, such as aerogel PHI HI]- However, for nematic 
elastomers, the network elasticity induces a long-range 
interaction between nematic directors. Furthermore, it 
is not a priori clear whether the effects of network het- 
erogeneities can be properly modeled as quenched ran- 
dom fields. Hence, the original Imry-Ma type of argu- 
ment does not necessarily apply. Nevertheless, absent 
any convincing alternative, this traditional random-field 
approach has been adopted by various authors [SHT3]. 

Simulation work has been done by Uchida [lOl [11] 
on a simple, two-dimensional, random-field model of ne- 
matic elastomers. He found that IGNEs always exhibit 
the PD state at a sufficiently low temperature. On the 
other hand, Terentjev et al. [9l [131 [IH have used Gaus- 
sian variational methods to study a similar model, and 
they concluded that replica symmetry is broken, indica- 
tive of glassy low-temperature properties in PD IGNEs. 
Xing and Radzihovsky [22j [23] start from a disorder- 
free model of a MD nematic elastomer, turn on an in- 
finitesimal quenched random field and study the con- 
sequent stability of the long-range nematic order, via 
renormalization-group techniques. They found that ne- 
matic order persists even below three spatial dimensions. 
The random-field idea is also implicit in the more re- 
cent works by Biggins et al. [T71 [TH], which employs a 
quasi-convexification method to construct textured, soft 
deformations of IGNEs, and thus allows these authors to 
obtain sharp bounds on the stress-strain curve. 

The aim of the present paper is to explore the nature 
of quenched disorder in IGNEs, and its impact on local 
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FIG. 1. (Schematic) snapshots of nematogen orientations at 
a particular instant (blue unshaded), and at a much earlier 
instant (shaded), (a) A conventional liquid crystal in the 
isotropic state. Such systems do not memorize the local pat- 
tern of nematogen alignment indefinitely. There is no correla- 
tion between the orientations of blue and shaded nematogens 
that are depicted near one another, (b) A liquid cystalline 
ealstomer in the macroscopically isotropic state. Such sys- 
tems do memorize the local pattern of nematogen alignment 
indefinitely. The orientations of blue and shaded nematogens 
depicted near one another are correlated. 



nematic order, within the framework afforded by vulcan- 
ization theory pm27j . We find that the most relevant 
type of disorder does indeed resemble a quenched ran- 
dom field, but with an important difference which we 
shall discuss in detail below. By using the replica Gaus- 
sian variational approach, we calculate the thermal as 
well as sample-to-sample fluctuations of the local nematic 
order. Our principal findings are that the PD state is al- 
ways stable, even at sufficiently low temperature and, 
furthermore, that for sufficiently strong disorder, the lo- 
cal frozen nematic order has an oscillatory character (i.e., 
anti-alignment) at short distances. 



II. ORDER PARAMETERS, FREE ENERGIES 
AND CORRELATORS 



The physics of randomly crosslinked systems depends 
on both the state under which the system is crosslinked, 
and the state under which it is measured, ft is therefore 
of crucial importance to distinguish between these two 
states. We call the former the preparation state and the 
latter the measurement state. The local nematic order 
Q° in the preparation state and its counterpart Q in the 
measurement state are, in principle, distinct, and this 
distinction is naturally captured by the formalism of vul- 
canization theory. For IGNEs crosslinked at high tem- 
perature, Q*^ is negligibly small. In contrast, although 
the macroscopic average of Q vanishes in IGNEs, its lo- 
cal value can be substantial; this would reflect random 
local frozen nematic alignment. It is also important to 
note the fundamental asymmetry in the causal ordering 
of Q" and Q: whilst Q'^ can influence Q, the converse 



is forbidden. In the framework of vulcanization theory, 
this causality is naturally enforced via the replica trick, 
in which the measurement-state order parameter is repli- 
cated n times, {Q";a = l,...,n}, and the number of 
copies n is later taken to zero. 

As IGNEs are qucnch-disordercd systems, we need 
to distinguish between thermal averaging (i.e., averag- 
ing over thermal fluctuations), denoted (•••), and disor- 
der averaging (i.e., averaging over the realizations of the 
crosslinking) , denoted [•••]. Assuming that the system is 
self-averaging, it follows that every experimentally mea- 
sured quantity has already been disorder-averaged. Var- 
ious disorder- and thermal-averaged quantities involving 
the local nematic order in the measurement state Q can 
be constructed. To lowest order in Q, one can construct 
the quantity [{Q{r))]; however, this quantity trivially 
vanishes, owing to the macroscopic isotropy of the PD 
state. At second order in Q, one can form the following 
two, distinct, averaged quantities: (i) the thermal corre- 
lator 

C^ir - r') Tr [( (Q(r) - (Q(r))) (Q(r') - (Q(r'))))] , 

(la) 

and (ii) the glassy correlator 

C'^(r-r') :=Tr[(Q(r))(Q(r'))]. (lb) 

Although we are focusing on the Cartesian scalar aspects 
of the correlators, it is straightforward to reconstruct the 
full structure of the corresponding fourth-rank Cartesian 
tensors, by appending suitable isotropic tensor factors 
constructed from Kronecker deltas. 

The glassy correlator automatically vanishes for 
isotropic liquids. For IGNEs, however, it is nonzero and 
encodes the correlation in the randomly frozen nematic 
pattern. The difference between these two states is qual- 
itatively indicated in Fig. [l] The value if the glassy 
correlator at the origin, C'^(r)|.p^Q, is the nematic ana- 
logue of the Edwards- Anderson order parameter for spin 
glasses pS], and measures the magnitude of the local 
frozen nematic order. We shall call this the glassy or- 
der parameter. On the other hand, the thermal correla- 
tor can be used to detect the existence of a continuous 
phase transition. If such a transition exists, the thermal 
correlator would diverge at the critical temperature. If a 
discontinuous transition occurs, however, the correlator 
would diverge at the spinodal point. 

To address the interplay between nematic ordering and 
network heterogeneity in randomly crosslinked systems, 
we shall also need the vulcanization order parameter 
57(r*', r^, . . . , r"), which is a scalar function of 1 + n 
replicated three-dimensional position vectors {r°']a — 
0, 1, . . . ,n} (see, e.g., Refs. [5iH?7| ). Up to an additive 
constant, it gives the joint probability density function 
(pdf) for the results of the following 1 -\- n indepen- 
dent measurements: At a time just prior to crosslink- 
ing, a monomer is found at r", whilst in n independent 
measurements made long after crosslinking, the same 
monomer is found at r^,...,r". This order parame- 
ter detects and characterizes the positional localization 
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n{f) = G d^z drpir) 



of monomers that emerges at the gelation/vulcanization 
transition. In the hquid phase, aU monomers are delocal- 
ized, so the joint pdf is a trivial constant, and hence the 
order parameter vanishes identically. In the gel/rubber 
phase, a fraction G of monomers are localized. Their 
positions in the 1 + n distinct measurements are then 
correlated with one another, and consequently the order 
parameter has a nontrivial, mean-field profile, given by 

27r/ F"' 

(2) 

where f := (r°, . . . , r") and V is the volume of the sys- 
tem. Note that the saddle point order parameter Cl is 
invariant under simultaneous translation of all replicated 
vectors f := (r", . . . , r"). The parameter r has the physi- 
cal meaning of the local inverse square localization length 
for the monomers, and it fluctuates from monomer to 
monomer with pdf p{t). In this work, we shall make 
the simplifying assumption that the localization length 
has a sharp value rather than being distributed, i.e., 
p{t) = 5{t — This does not change our essential 

results. 

It is rather significant that the heterogeneous nature 
of elastomers is captured already at the mean field level 
of vulcanization theory. At long lengthscales, the most 
important coupling between the vulcanization order pa- 
rameter n at the saddle point and the nematic order 
parameters Q" reads 

/n 
df ^vrV^"n(f)vfvfn(f)Q«(r")Qf,(r'^)(3) 

where df indicates multiple integration over all replicated 
positions r" , and a summation is implied over the Carte- 
sian coordinates .... In the above summation, the 
replica-diagonal terms and the replica-off-diagonal terms 
contribute qualitatively distinctly from one another. In 
particular, the replica-diagonal terms can be omitted 
as they simply renormalize the quadratic term in the 
Landau-de Gennes free energy We proceed to in- 

tegrate out the coordinates r"' (with 7 ^ ct,f3). Hence, 
Eq. ([3| is transformed into the following, essentially non- 
local, quadratic term: 

/ d\' uj{r ~ r') J2 Tr Q"(r-) Q'^{r'), (4a) 



,li = 



(4b) 



where G is the gel fraction, and w is a kernel that rep- 
resents the disordering effects of the random network on 
the nematic order. The characteristic lengthscale beyond 
which uj{r) is suppressed is S^l- This suppression reflects 
the liquid nature of the network at short lengthscales but 
not at lengthscales longer than S^l. 

For IGNEs crosslinked at a very high temperature, the 
local nematic order Q*' in the preparation state is negligi- 
bly small and can therefore be ignored. The full free en- 
ergy of IGNEs is a combination of the Landau-de Gennes 



free energy for local nematic order in the measurement 
state and the coupling term, Eqs. 

H[Q]^J d\ |^TrQ"Q" + :|TrV,Q"V,Q" 



^Tr(Q")'V|(TrQ"Q")^ 

h 
2 



(5) 



/dVdVa;(r-r')Tr Q"(r)Q'^(r') 



,/J=r 



Here, T,. is the reduced temperature in the measurement 
state. In addition, K is proportional to the Frank con- 
stant in the one-constant approximation. It is important 
to note that any standard replica treatment of a phe- 
nomenological theory in which the nematic order is cou- 
pled linearly to a quenched Gaussian random field with 
correlator w(r) would necessarily lead, in addition to the 
replica off-diagonal term in Eq. ([5|, to a nonvanishing 
replica-diagonal term having precisely the same kernel as 
its off-diagonal counterpart. This absence of the replica- 
diagonal term is a salient feature of vulcanization theory. 
It demonstrates that the quenched disorders in elastomers 
can not be simply modeled as random field. As we shall 
show below, it is also responsible for the oscillatory spa- 
tial structure for IGNEs in the strong-disorder regime. 
Standard random-field models would not yield such os- 
cillatory structure. 

Focusing on the quadratic parts of the effective Hamil- 
tonian ([5|, we can obtain the bare thermal {Cq) and 
glassy (C^) correlators [see Eqs. ([l])] of the nematic order 
parameter in wave- vector space: 



Ri 



t + K^ + Ae-^y^' 

RiAer^"/^ 
(t + K^ + Ae-'^Vz) 



(6a) 
(6b) 



Here, Ri 



._ f2 



^i/K, whilst t 



Tra/K, A 



32TT^hG^ / K^j^ > and k := pS^i^ are the dimensionless 
rescaled measurement-temperature, disorder strength, 
and wavevector. The appearance of A in the denomi- 
nators, in combination with the scale- dependent factor 
exp(— K^/2), is a significant result of the present ap- 
proach. It follows directly from the absence of the replica- 
diagonal contributions to the term Q, and therefore 
cannot result from a standard random-field approach. 
Originating in the random network, and proportional 
to the shear modulus of the network [5D], it leads to 
a scale-dependent, downward renormalization of the re- 
duced temperature. 

The bare correlators exhibit instability at some neg- 
ative reduced temperature t < 0. For A < 2, i.e., at 
weaker disorder, the denominator factor t + K^ + Ae~'^ 
has a minimum at k = 0. Consequently, both of the cor- 
relators given in Eqs. Q diverge, at k = 0, at a criti- 
cal (reduced, rescaled) temperature — A. This suggests a 
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continuous transition to a uniformly ordered (i.e., macro- 
scopically anisotropic) nematic state. By contrast, for 
A > 2, i.e., at stronger disorder, the denominator fac- 
tors have a minimum at = 21n(A/2); therefore, both 
bare correlators exhibit nonzero-wavelength divergences 
at a (reduced, rescaled) temperature — 21n(eA/2). This 
suggests a continuous phase transition to a state having 
periodic modulations. As we shall see shortly, however, 
neither of these putative transitions actually occurs, both 
of them being precluded by the quenched randomness. 



III. GAUSSIAN VARIATIONAL APPROACH 
TO THERMAL AND GLASSY CORRELATORS 

To investigate the fate of the system at temperatures 
below these putative instabilities, we employ the Gaus- 
sian variational method (see, e.g., Ref. [31]) for systems 
having quenched randomness 32^35]. We choose the cor- 
responding quadratic trial Hamiltonian Hq to be 



By minimizing i^var with respect to and C'^ for a given 
value of Q, and defining a renormalized reduced temper- 
ature t]i via 



= t + 7wRi J ( C 



1 



TrQ^ 



(11) 



we obtain for the renormalized thermal and glassy corre- 
lators the results 



Ri 



RiAe-'^'^^ 



{tR + K^ + Ae-'^y^y 



(12a) 
(12b) 



Note that Eqs. ( |12| are structurally identical to their bare 
counterparts, Eqs. (|6|, the only difference being the re- 
placement of the bare reduced temperature t by its renor- 
malized counterpart tR. 



The relation (11) between tR and t determines the ex- 



istence, or lack thereof, of a transition out of the macro- 
scopically isotropic state. Let us first discuss the high 
(measurement) temperature regime, in which Q = and 
the correlators are small. In this case, the second term on 



Tr(Q"(pi) Q(pi))(Q (P2) Q(P2)), (7) the RHS of Eq. 0, proportional to w, can be neglected 



which generically features a nonvanishing mean nematic 
order parameter Q together with a kernel F"^ Here, 
the notation denotes / d^p/ {2t:Y. Assuming that 

replica permutation symmetry remains intact I37j . we 
may parametrize the replica-space inverse of T"" as 



Thus, tR is essentially the bare (reduced, rescaled) tem- 
perature t, and the renormalized correlators ( 12 ) reduce 
to the bare correlators ([6]), as expected. 

But what happens at lower temperatures? As t is re- 
duced, the correlators become larger, and the corrections 
due to fluctuations, present in Eq. (11), increase. The 



r^^)"'^(Pi P2) '— '^(12) (C"^(P]^)(5"'^ +C'^{pi)') (8) essential question is then the following: Do the denonii- 



where 5(12 • • • m) denotes {2itY8{pt^ +P2 + ' ' ' +Pm)' 

(jp) and C^{p) are the renormalized thermal and glassy 
correlators, in which fluctuations are approximately ac- 
counted for. Both of them, together with Q, are to be 
determined self-consistently, by minimizing the resulting 
variational free energy 

/ n 
l[VQ'-e-"°, (9) 
a=l 

where {■■■)ho denotes averaging with respect to the 
Boltzmann weight e^^" . The explicit form of Fvar is 



nators in Eqs. (12) vanish at some wave-vector p for low 
enough temperature? If this were indeed to happen, it 
would signify a continuous transition to some macroscop- 
ically symmetry-broken state. Now, the low-temperature 
physics of our model depends sensitively on the value of 



A. For A < 2, the denominator factor in Eqs. (12) [viz., 

tR + K,"^ + Ae"" has a minimum at k = 0. Assum- 
ing that a continuous transition occurs at some finite 
temperature tc, this denominator would have to vanish 
for K = and tR + A = 0. For three-dimensional sys- 
tems, however, the integral of the glassy correlator /^C^ 
in Eq. (11) has infrared divergence, and thus tR{tc) di- 



5nV 



2\nT 




Ae-'^'/2)(C^ + TrQV5y) 



5V. 



Pi 



Tr / (C + C- 



P 



P2 



5(123)TrQ^ 



verges, implying that the denominator would diverge, 
and this contradicts the assumption that this denomi- 
nator vanishes. We therefore conclude that in three di- 
mensions and for A < 2 there is no continuous transi- 
tion to any state having long-range nematic order. For 
A > 2, the factor tR-\- + Ae~'^ has a minimum at 
— 2 ln(A/2). In this case, to have a continuous transi- 
tion (towards a macroscopically periodically modulated 
state) would require the vanishing of the denominator, 
now at = 21n(A/2). However, once again, the integral 



PlP2'P3 



in Eq. ( 11 1 is divergent, and therefore no continuous tran- 



5(1234) (TrQ2 



(10) 



sition occurs. Hence, owing to network heterogeneities, 
the macroscopically isotropic phase is always locally sta- 
ble, for any positive A. This result agrees with experi- 
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FIG. 2. (a) Glassy order parameter C'^(r)|^^Q (measured in 
unit Ri) and (b) glassy correlation length ^G,d (in unit (^l) as 
a function of t for three cases of weak disorder: A = 0.2 (red 
solid), A = 0.7 (green dashed) and A = 1.5 (blue dot-dashed), 
where we have set wRi = OA-k. 



mental observations, e.g. by the Urayama group |16j . and 
the results of simulations performed by Uchida [TOl [H] . 



IV. CORRELATOR STRUCTURE AND 
LOW-TEMPERATURE PROPERTIES 

We now unfold the physics encoded in the correlators 



Eqs. (12 1. First, we look at the glassy order para mete r 
C'^(r)]7lo' which is given by J^C^{p). From Eq. (12b|, 
we see that below t = it increases rapidly with decreas- 
ing temperature, becoming, at low temperatures, asymp- 
totically linear in t, behaving as i4r|t|/7w^|. We have 
computed this order parameter: the results are shown 
for three cases of weak disorder in Fig. [2^ and one case 
of strong disorder in Fig. |4)d. 

As the fluctuation-response theorem indicates, the 
static susceptibility x{p) for local nematic order is related 
to the thermal correlator via x{p) = T~^C^ (p). The 
peak value of , which for weak disorder (i.e., A < 2) 
occurs at = and for strong disorder (i.e., A > 2) 
occurs at = 21n(A/2), increases rapidly as the tem- 
perature is decreased, as we see from Eq. (12a). Due to 
the coupling between local nematic alignment and elastic 
strain, we thus expect there to be a corresponding soft 
elasticity. We speculate that this mechanism lies behind 
the supersoft elasticity of PD IGNEs (cf. Refs. [TTlfTS] ). 
A detailed analysis of elasticity of PD IGNEs is however 
not possible using the current simple model which does 
not involve the elastic degrees of freedom. We shall pur- 
sue it in a separate publication. 

In the weak-disorder regime (i.e. A < 2), both corre- 
lators have a finite peak at zero wave-vector, indicating 
spatial decay (but without oscillation). The characteris- 
tic lengthscales £^T,d and ^G,d over which the correlators 
decay can be obtained by examining their small wave- 
vector behaviours: 



(2- 



2{tR + A) 



^G,d 



'^^T,d 



lf2 



(13) 



Since t^; -I- A approaches zero towards low temperature, 
both lengthscales ^G/T,d grow rapidly (see Fig.jljD). 




FIG. 3. (a) Real-space optical microscopic image of polydo- 
main structure in IGNEs, obtained by the Urayama group. 
Note the occurrence of a short-range 'checkerboard' pattern 
indicative of oscillatory structure in the local nematic align- 
ment, most evident in the upper-right region of the panel, 
(b) Noise-filtered (Log amplitude-squared) Fourier transfor- 
mation of the optical image in panel (a). Note the occur- 
rence of a central high-intensity (bright) peak accompanied 
by a high-intensity ring, which correponds to the checker- 
board pattern in real space. The shaded rectangular grid is 
due to the denoising filter. 



In the strong-disorder regime (i.e., A > 2), the ther- 
mal correlator and the glassy correlator C*^ are each 
peaked on a spherical shell of wave- vectors: C^{p) peaks 



for = 21n(A/2)/^2^ and peaks for \p\' 
with ^G,o obeying 

ifl + 4+(a/eG,o) 



_^e-{U/ia,.r/2 (14) 



We term ^g,o the glassy oscillation length. The peak- 
ing of C'^{p) and C*^(p) is associated with oscillatory 
behaviour (i.e., anti-alignment) of the local nematic or- 
der, which is superposed on their overall decay in real 




FIG. 4. Glassy correlator C'^ (measured in unit Ri) at A = 
10 as a function of reduced temperature t and (a) reduced 
wave- vector k; (b) separation distance r (in unit ^l), with 
wRi = 0.47r. Note the progression in (a) of the peak location 
from zero wave- vector to nonzero wave- vector with decreasing 
t. The peak at nonzero wave-vector is associated with the 
spatial oscillations observed in (b) for the glassy correlator in 
real-space. Note also in (b) that the glassy order parameter 
C'^{r)\r^o rapidly increases as the temperature is reduced. 
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FIG. 5. (a) Crossover diagram for the glassy and thermal 
correlators, indicating the three qualitatively distinct regimes 
of behavior. Above the blue solid line, both correlators oscil- 
late and decay as a function of separation. Between the blue 
solid and red dashed lines, both correlators decay but only the 
thermal one also oscillates. Below the red dashed line, both 
correlators decay but neither oscillates, (b) Oscillation wave- 
length (,c,o (measured in unit ^i) as a function of reduced 
temperature t for A = 10. As t decreases, ^g,o eventually 
saturates at a nonzero value of order For both panels, we 
have set wRi = OAn. 



space. This echos the polydomain structure with well de- 
fined length scale observed by the Urayama group [IB] . 
Indeed, as shown in Fig. [3j Fourier transform of their 
real space image of the PD structure, after denoising 
and deblurring, shows clearly a ring of local maximum 
at the scale of about half micron. In Fig. |4j we show 
examples of how the shell radius for the glassy correla- 
tor progresses from the origin to a nonzero wave-vector 
magnitude and, furthermore, how the oscillatory decay 
develops in real space, as the measurement temperature 
is decreased. The crossover, from non-oscillatory to os- 
cillatory decay as the disorder strength is increased, is 
shown in Fig. [5^. It can be deduced from Eq. (14 1 that, 
as the temperature is decreased, ^g,o decreases and ul- 
timately saturates at the temperature-independent ther- 
mal oscillation length ^t,o (= Ci/\/21n(A/2)). Provided 
that ^G,d exceeds £,g,o, the length „ will reflect the size 
of the locally aligned nematic domains; therefore, the de- 
crease of ^G.o (shown in Fig.[5j3) will reflect the shrinking 
of the domain size with decreasing temperature. 



V. CONCLUSION AND OUTLOOK 

Via a strategy rooted in vulcanization theory, we have 
developed a unified approach to the statistical physics 



of nematic elastomers, in which nematic ordering — both 
in the preparation and the measurement states — and net- 
work heterogeneity are naturally incorporated. It is note- 
worthy that the quenched disorder naturally emerges 
from this approach differs in essential ways from the tra- 
ditional quenched random field approach. This differ- 
ence evinces a strength of vulcanization theory. Using 
a variational method, we show that the macroscopically 
isotropic state of isotropic-genesis nematic elastomers is 
always thermodynamically stable and, furthermore, that 
for sufficiently high disorder strengths, short-range oscil- 
latory structure of the local nematic alignment arises. 

Apart from its relevance to the specific subject of liquid 
crystalline elastomers, the present work brings to light a 
more general issue, viz., that the concept of a quenched 
random field should be broadened to incorporate not only 
the traditional, 'frozen' type, which does not fluctuate 
thermally, but also the type necessary for understand- 
ing media such as liquid crystalline elastomers, in which 
the frozen nature of the random field is present only at 
longer lengthscales, fading out as the lengthscale pro- 
gresses through a characteristic localization length, owing 
to the thermal position fluctuations of the network's con- 
stituents. The framework elucidated in the present work 
can be extended, with suitable modifications, to explore 
the statistical physics of other randomly crosslinked sys- 
tems, such as smectic elastomers and various biological 
materials. 



ACKNOWLEDGMENTS 



We thank Kenji Urayama for informative discus- 
sions on experimental results and Xiaoqun Zhang for 
valuable assistance with image processing. This work 
was supported by the National Science Foundation via 
grants DMR06-05860 and DMR09-06780, the Institute 
for Condensed Matter Theory and the Frederick Seitz 
Materials Research Laboratory at the University of Illi- 
nois at Urbana-Champaign, the Institute for Complex 
Adaptive Matter, and Shanghai Jiao Tong University. 



[1] Warner, M. & Terentjev, E. M. Liquid Crystal Elas- 
tomers. Oxford University Press (2003). 

[2] Olmsted, P. D. Rotational invariance and Goldstone 
modes in nematic elastomers and gels. J. Phys. II 
(France) 4, 2215-2230 (1994). 

[3] Golubovic, L. & Lubensky, T. C., Nonlinear elasticity of 
amorphous solids Phys. Rev. Lett. 63, 1082-1085 (1989). 

[4] Verwey, G. & Warner, M. Soft rubber elasticity. Macro- 



mol. 28, 4303-4306 (1995). 
[5] Ye, F., Mukhopadhyay, R., StenuU, O. & Lubensky, T. C. 
Semisoft nematic elastomers and nematics in crossed 
electric and magnetic fields Phys. Rev. Lett. 98, 147801 
(2007). 

[6] Deam, R. T. & Edwards, S. F. The theory of rubber 

elasticity. Phil. Trans. R. Soc. A 280, 317-353 (1976). 
[7] De Gennes, P. G. Scaling Concepts m Polymer Physics. 



7 



Cornell University Press, Ithaca (1979). 
[8] Clarke, S. M., Terentjev, E. M., Kundler, I. & Finkel- 
mann, H. Texture evolution during the polydomain- 
monodomain transition in nematic elastomers. Macro- 
mol. 31, 4862-4872 (1998). 
[9] Fridrikh, S. V. & Terentjev, E. M. Polydomain- 
monodomain transition in nematic elastomers. Phys. 
Rev. E60, 1847-1857 (1999). 

[10] Uchida, N. Elastic Effects in Disordered Nematic Net- 
works Phys. Rev. E 60, R13-R16 (1999). 

[11] Uchida, N. Soft and nonsoft structural transitions in 
disordered nematic networks. Phys. Rev. E 62, 5119- 
5136 (2000). 

[12] Selinger, J. V. & Ratna, B. R. Isotropic-nematic transi- 
tion in liquid-crystalline elastomers: lattice model with 
quenched disorder. Phys. Rev. E 70, 041707 (2004). 

[13] Petridis, L. & Terentjev, E. Nematic-isotropic transition 
with quenched disorder. Phys. Rev. E 74, 051707 (2006). 

[14] Petridis, L. & Terentjev, E. Quenched disorder and spin- 
glass correlations in XY nematics. J. Phys. A: Math. 
Gen. 39, 9693-9708 (2006). 

[15] Feio, G., Figueirinhas, J. L., Tajbakhsh, A. R. & Teren- 
tjev, E. M. Critical fluctuations and random-anisotropy 
glass transition in nematic elastomers. Phys. Rev. B 78, 
020201 (2008). 

[16] Urayama, K., Kohmon, E., Kojima, M. & Takigawa, T. 
Polydomain-monodomain transition of randomly disor- 
dered nematic elastomers with different cross-linking his- 
tories. Macromol. 42, 4084-4089 (2009). 

[17] Biggins, J. S., Warner, M. & Bhattacharya, K. Supersoft 
elasticity in polydomain nematic elastomers. Phys. Rev. 
Lett. 103, 037802 (2009). 

[18] Bhattacharya, K., Biggins, J. S. & Warner, M. 
Elasticity of polydomain liquid crystal elastomers. 
\arXw:0911. 3513^ 1, 2009. 

[19] Imry, Y. & Ma, S.-K. Random-field instability of the 
ordered state of continuous symmetry. Phys. Rev. Lett. 
35, 1399-1401 (1975). 

[20] Radzihovsky, L. & Toner, J. Smectic liquid crystals in 
random environments. Phys. Rev. B 60, 206-257 (1999). 

[21] Feldman, D. E. Quasi-long-range order in nematics con- 
fined in random porous media. Phys. Rev. Lett. 84, 4886- 
4889 (2000). 

[22] Xing, X. & Radzihovsky, L. Universal elasticity and fluc- 
tuations of nematic gels. Phys. Rev. Lett. 90, 168301 
(2003). 

[23] Xing, X. & Radzihovsky, L. Nonlinear elasticity, fluc- 
tuations and heterogeneity of nematic elastomers. Ann. 
Phys. 323, 105-203 (2008). 



[24] Goldbart, P. M., Castillo, H. E. & Zippelius, A. Ran- 
domly crosslinked macromolecular systems: vulcaniza- 
tion transition to and properties of the amorphous solid 
state. Adv. Phys. 45, 393-468 (1996). 

[25] Peng, W., Castillo, H. E., Goldbart, P. M. & Zippelius, 
A. Universality and its origins at the amorphous solidi- 
fication transition. Phys. Rev. B 57, 839-847 (1998). 

[26] Broderix, K., Weigt, M. & Zippelius, A. Towards finite- 
dimensional gelation. Eur. Phys. J. B 29, 441-455 
(2002). DOI: 10.1140/epjb/e2002-00325-4 

[27] Xing, X., Pfahl, S., Mukhopadhyay, S., Goldbart, P. M. 
& Zippelius, A. Nematic elastomers: from a microscopic 
model to macroscopic elasticity theory. Phys. Rev. E 77, 
051802 (2008). 

[28] Edwards, S. F. & Anderson, P. W. Theory of spin glasses. 
J. Phys. F: Metal Phys. 5, 965 (1975). 

[29] For a = /3, we can integrate out all cc^ for y a. 
The V"V"nV"V"n part in Eq. ^ then reduces to a 
function of x" only. Due to the invariance of the saddle 
point order parameter Q under simultaneous translation 
of f := (r", . . . , r"), however, this function must be a 
trivial constant. 

[30] At the mean field level, the shear modulus of a network 
is proportional to G^^^'^, and so is the disorder strength 
A. 

[31] Feynman, R. P. Statistical Mechanics: A Set of Lectures. 

Westview Press, 2nd Edn. (1998). 
[32] Shakhnovich, E. I. & Gutin, A. M. Formation of mi- 

crodomains in a quenched disordered heteropolymer. J. 

Phys. (France) 50, 1843-1850 (1989). 
[33] Mezard, M. & Parisi, G. Interfaces in a random medium 

and replica symmetry breaking. J. Phys. A: Math. Gen. 

23, L1229-L1234 (1990). 
[34] Mezard, M. & Parisi, G. Replica field theory for random 

manifolds. J. Phys. I (France) 1, 809-836 (1991). 
[35] Stepanow, S., Dobryrnn, A. V., Vilgis, T. A., & Binder, 

K. Copolymer melts in disordered media. J. Phys. I 

(France) 6, 837-857 (1996). 
[36] Strictly speaking, the structure of the kernel F"'' should 

allow for anisotropy for the case of Q 7^ 0. Our main 

analysis, however, only concerns the isotropic state, for 

which Q = 0. 

[37] According to the standard understanding, a replica sym- 
metry breaking at low temperature implies that the sys- 
tem freezes into one of many nearly degenerated states. 
Such a system should exhibit irreversible glassy proper- 
ties such as hysteresis in stress-strain curve. These prop- 
erties have not been observed in polydomain nematic 
elastomers. 



